Phenomenological model for the Drell-Yan process: Reexamined 
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Drell-Yan pair production is investigated. We reexamine a model where the quark momentum fraction is 
defined as the ratio of the corresponding light cone components of the quark and parent nucleon in a naive 
parton-model approach. It is shown that results differ from the standard parton model. This is due to unphysical 
solutions for the momentum fractions within the naive approach which are not present in the standard parton 
model. In a calculation employing full quark kinematics, i.e. including primordial quark transverse momentum, 
these solutions also appear. A prescription is given to handle these solutions in order to avoid incorrect results. 
The impact of these solutions in the full kinematical approach is demonstrated and compared to the modified 
result. 

PACS numbers: 12.38.Qk, 12.38.Cy, 13.85.Qk 



I. INTRODUCTION 

The Drell-Yan (DY) process [1] was first described in the 
1970s and provides an important tool to access the distribu- 
tion of partons inside the nucleon. While a lot of informa- 
tion can be gained from deep inelastic scattering [2], mea- 
surements of Drell-Yan events give complementary insights, 
especially about sea-quark distributions [3]. This has sparked 
many studies of this process 10-01 which are generally in- 
spired by perturbative QCD (pQCD). A lot of experimental 
effort is being devoted to measurements of the DY process: In 
antiproton-proton collisions at PANDA (FAIR) 10] and PAX 
Jl, in protonjroton collisions at RHIC J-PARC HH- 

LLJ, IHEP Cj and JINR 0] and in pion-nucleon collisions 
at COMPASS |HE[i3l- An overview of the experimental sit- 
uation can be found in lfl8ll . PANDA, for example, will allow 
measurements at hadron cm. energies of a few GeV, where 
non-perturbative effects are expected to become more impor- 
tant. This highlights the need to model these effects in a phe- 
nomenological picture. 

In addition the standard pQCD leading order (i.e. parton 
model) description does not fully describe the interesting ob- 
servables. Invariant mass (M) spectra of the DY pair can 
only be accounted for by including an additional K-factor 
and transverse momentum (pj) spectra are not accessible at 
all IU9I1 . The latter can be partly cured by folding in a phe- 
nomenological Gaussian distribution for the transverse mo- 
mentum, the width of which has to be fitted to data. However 
the absolute size of the cross sections is still underestimated 
1I20I1 . Next-to-leading-order (NLO) calculations improve the 
description in some aspects, but also bring about additional 
problems. The calculated invariant mass spectra come closer 
to the data and the pj spectra are comparable to data in the 
region p T ~ M (hJ, but not for p T -» H. In fact the 
pr-spectrum is divergent for pj — > in any fixed order of 
the strong coupling a s , due to large logarithmic corrections 
In (M/ pj). These stem from soft gluon exchange and it is pos- 
sible to remove these divergencies by an all-order resumma- 
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tion. However since pj is no longer a hard scale at pj — > ad- 
ditional non-pertubative (i.e. experimental) input is needed in 
these (and all other pQCD) approaches to describe the region 
of very small pr 1121142311 . Note here that the parton model 
(i.e. leading order) description is still a very useful starting 
point, e.g. for studying spin asymmetries in DY, since there 
NLO corrections appear to be rather small I24U26I1 . 

A phenomenological model that incorporates full trans- 
verse momentum dependent quark kinematics and which in 
addition allows for mass distributions of quarks was proposed 
to resolve these problems l27l - l29ll . The idea stems from the 
fact that in the usual collinear approach the parton momenta 
are confined to the beam direction, thus only one momen- 
tum component is different from zero. The other components, 
namely the transverse momentum and the mass of the parton, 
do not enter into the calculation of the partonic subprocess 
cross section. Since at finite energies these components might 
influence the cross section to some extent it is worthwhile to 
examine this influence in detail. However it turns out that in 
these works important physical constraints were not consid- 
ered and thus incorrect results were obtained. In the current 
paper we examine these constraints in detail and present a pre- 
scription to properly account for them. Finally we compare 
the results of the treatment in 1271 - 12911 with our corrected re- 
sults. Since the mentioned problems already appear for the 
case without mass distributions for partons we restrict our- 
selves here to massless partons. 

This paper is organized as follows: in Sec. [TT] we com- 
pare the standard collinear parton model description for DY 
with an approach that defines the parton momentum fraction 
x via light cone components. The latter approach will be a 
demonstration of the problems that appear in the calculation 
with the full kinematics. Section|III]contains two calculations 
in the full kinematical scheme, i.e. taking into account the full 
transverse momentum dependence of the partonic subprocess. 
The approach of I27I - I29I1 is discussed in detail in Sec. IIIIB II 
and the technical details are given in the Appendix. Section 
IIII B 2l then contains our calculation which respects the physi- 
cal constraints laid out in Sec. MB I The numerical results are 
presented in Sec. |IV]where we compare the two calculations 
of Sec. [Ill] quantitatively. Finally we present our conclusions 
inSec.|V] 

In the following we present the conventions and notations 
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used throughout this paper: It will turn out to be useful to 
write four-momenta using light-cone coordinates. We employ 
the following convention for general four- vectors a and b 



a = ciq + a z , 
cT — Go — a z , 

a± = [a x , fl v ) , 



a - a a 



• b — — {a + b + a b + — 2 a ± ■ b±j . 



(1) 

(2) 

(3) 
(4) 

(5) 



We regard all particles as massless. We define the target nu- 
cleon to carry the four-momentum P\ and the beam nucleon 
to carry the four-momentum Pi (see fig. [TJ. In the hadron 
center-of-mass (cm.) frame we choose the z-axis as the beam 
line and the beam (target) nucleon moves in the positive (neg- 
ative) direction. Therefore the nucleon four-momenta read 



Pi = 



( 2 2 

2 2 



(6) 
(7) 



Note here that with a finite nucleon mass Pi and P2 would 
change. We have explicitly conducted the entire calculation 
with non-zero nucleon mass and convinced ourselves that it 
does not influence our arguments in Sees. [ID and [Til] Our 
results in Sec. [TVjwould receive only a small correction, since 
we are looking at cm. energies of y[S > 27 GeV. Thus, in 
this paper, we have put the nucleon mass to zero for the sake 
of simplicity and readability. 

We denote the four-momentum of the parton in nucleon 1 
(2) as p\ (pi)- The on-shell condition in light-cone coordi- 
nates then reads 



O-pl^PtPi-iPiJ 1 



The definition of the Feynman variable xf is [3 

1z 



X F 



(q-X 



(8) 



(9) 



For the virtual photon in fig.[T]the maximal q z is derived by re- 
quiring the invariant mass of the undetected remnants to van- 
ish and the photon to move collinearly to the nucleons: 



(/>! + P 2 ~ q) 2 = X 2 = 



(10) 



5+<? 2 -2V5 r A /^ + (^ ax = (11) 
S -q 2 



2yfS 



(qz)n 



(12) 



II. COLLINEAR APPROACH 




FIG. 1. DY production in a nucleon-nucleon collision; X[ and Xi 
denote the nucleon remnants. See main text for details. 



parton model with a naive approach which uses the light-cone 
component definition of the parton momentum fractions. It 
will turn out that in the latter case unphysical solutions appear 
that must be removed to be consistent with the standard parton 
model. 



A. Standard parton model 

The leading-order Drell-Yan total differential cross section 
in the standard parton model reads |Q]] 



dcr = J djt] J dx2^q 2 fi(xi,q 2 )f 1 (x 2 ,q 2 )dd-(xi,X2,q' 

(13) 

Here x\ and X2 are the momentum fractions carried by the 
annihilating partons inside the colliding nucleons: 



Pi 

P2 



■■XiPi 

x 2 P 2 



(14) 
(15) 



The sum runs over all quark flavors and antiflavors, q, denotes 
the electric charge of quark flavor i, the functions /; are parton 
distribution functions (PDFs) and dcr is the total differential 
cross section of the partonic subprocess, 



Ana 1 
"V" 



5(M 2 - q 2 ) 6 (A \pi +p 2 -q) A dM 2 . (16) 



Here q is the four- momentum of the virtual photon, P\,P2 are 
the four-momenta of the partons (cf. fig.Q]) and a * 1/137 is 
the fine-structure constant. 

Note that it becomes immediately clear from Eqs. ( TT4l > and 
([T5l l that the incoming partons move collinearly with the nu- 
cleons. According to Eq. (TT~6b no transverse momentum can 
be generated for the virtual photon (and thus for the DY pair) 



In this section we treat the interacting partons as collinear 
with their parent nucleons. We compare the standard textbook 
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in the leading-order process: The maximal information about the DY pair that can be 

gained from Eq. ( fT3l is double differential. A common choice 
Pi ± — P2 ± = (17) of variables is the squared invariant mass M 2 and Feynman's 

=> $ i4 \pi +P2~q) = ff((pl + p2? ~ q°) 6 (2 \q±) Xf ° f the viltUal P hot ° n: 

X ff((pi + Plf - q z ) ■ (18) 



d(3- 



dM 2 dx f 



I 



Ana 



d A q—5{M 2 -q 2 )6 m {p l + Pl -q)5\x F 



q : 



(ft)n 

r 



Ana 2 
9M 2 



5{M 2 ~( Pl + P2 f) d(x F 



(Pi + Pz)z 



(19) 



The two ^-functions connect x\ and x 2 with the chosen ob- 
servables 



M = 2p\pi = X\X%S 

VS (x 2 - Jti) 



xf = +- 



2(ft) n 



with fe) max = cf. Eq. 03- 

Solving for X\ and x 2 yields 



*U = 



max XF + y]((q Z )maxXF) 2 +M 2 



— (^fz)max -^F — ^xoll 

Vf ' 



(20) 
(21) 



(22) 



A '2 ± 



(ft) max 



(ft) max 

with the energy of the collinear DY-pair 



^coii = tJm 2 + ((q z ) max x F ) 2 . 



(23) 



(24) 



However the lower solutions are always negative. Only the 
upper solutions are in the integration range of Eq. ( TT3l and are 
physically meaningful. For the negative solutions the parton 
energies would be negative on account of Eqs. ( TBI and ( fT5l ). 
The hadronic cross section then reads: 



d<T 



= ^ dx 2 J] qlftx, , M 2 ) fi(x 2 , M 2 ) ^S(M 2 - (p, + p 2 ) 2 ) 5 - 

= f d*i f di 2 y^(i,,M 2 )/;(i2,M 2 )^ (gz)max ^(x l -x l+ )^(x 2 -x 2+ ) 

Jo Jo ' 9M S^q z )L x x 2 F+ M 2 

= V qHix^M^Mx^M 2 ) ^4 ■ (25) 



In this section we have presented the standard parton model 
solution for the leading-order DY cross section. The only 
quantities in this approach not determined by pQCD are the 
PDFs. These have to be obtained by fitting parametrizations to 
experimental data, mainly on deep inelastic scattering (DIS), 
but also on measurements of DY production itself IB ill . 

B. Naive parton model 

In this section we work out the complete collinear kinemat- 
ics using the definition of the parton momentum fraction as the 



ratio of light-cone components of the parton and the nucleon 
Hit]. We show that there exist other solutions for the parton 
momentum fractions x, which are neglected in the standard 
parton model right from the start. These other solutions will 
turn out to be unphysical and are derived at this point only 
to provide insight into difficulties arising from a transverse- 
momentum dependent calculation as discussed in Sec. IIIIB II 
The partons inside the nucleons carry some fraction of their 
parent hadron's longitudinal momentum. Labeling the par- 
ton momentum inside nucleon i with /?, we can define these 
fractions as ratios of plus or minus components of the par- 
tons and the corresponding components of the nucleon mo- 
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menta. In the Drell-Yan scaling limit (S — > oo and M 2 /S 
finite) P~ = P^ — VS~ become the large components while all 
other components vanish. Note here that with a finite nucleon 
mass the large components would be modified and the small 
components would be nonzero. This however poses no prob- 
lem for the following calculations, c.f. the discussion below 
Eqs. © and ©. We define 



Xl = 



Xz = 



1 

P l 

n~ 

Pi 



Pi 



Pj_ 



(26) 
(27) 
(28) 
(29) 



Note that Eqs. (|26j> and (|28j are standard definitions JH]. The 
tilde quantities in Eqs. (|27T i and d29l are introduced for later 
convenience. The kinematical constraints for these fractions 
are the on shell conditions 



pi 



p\p\ 
p\pi 



: 
: 



x\X\- . 
xixi- . 



together with 



M 2 = (pi +pi) 2 = 2pxp 2 = PiP 2 +PiPi 
— (x\Xi + xix 2 ) S 



(30) 
(31) 



(32) 



and 



Xf 



(Pi +Pi) z 



l 



2(<?z)ma 

_Js_ 

2(<7.-)ma 



■(Pl-Pl +Pl~Pl) 

■ (x~\ -x\+xz - x 2 ) . 



(33) 



We will show now that the constraints in Eqs. <f30b-(f33l> can 
be fulfilled by two different sets of momentum fractions x ; , x t . 
Equation ( 1301 implies x\ — or x\ = 0. If 



xi = 

Eg. Qll M 2 

=> xi + + xi 



X\Xl 



Eq. HU 



Eq. {33} 





(x 2 - Xl) 



V5 



2(? z )n 



(34) 

(35) 
(36) 
(37) 

(38) 



This is just the standard parton-model solution, Eqs. ( 1201211 ). 
as described in Sec. Ill Al However there exists another solu- 



tion, namely for x\ = 0: 



xi =0 

Eq. ED M 2 



X\Xl 



X~l+§ + Xl 



Eq. E} 



xi =0 



Eq.(33} 

> X/r = (Xl - X2); 



2(^)n 



(39) 

(40) 
(41) 
(42) 

(43) 



Kinematically this second solution represents the (strange) 
case where each parton moves into the opposite direction of 
its respective parent nucleon. One can see this in the following 
example, where we choose x$ = 0. Then we have 



xi = x 2 



M 



Pi = 



M 



\{P\-P,)-\^M 2 



and analogously 



Pi 



M 
~~2 



(44) 
(45) 

(46) 



Since nucleon 1 (2) moves into negative (positive) z-direction, 
cf. Eqs. ((6]) and (0, the partons here move exactly opposite. 
The parton momentum fractions x, (not x, !) entering the PDFs 
in Sec. Ill Al however are those of partons that move into the 
same direction as their parent nucleon. The second solution is 
thus physically not meaningful and it is discarded right away 
in the standard parton model approach. 

The essential difference between the standard and the naive 
parton model is the following: In the (collinear) standard par- 
ton model all components of p, are fixed at once by pi = XiP u 
This automatically implies x,- = 0. Such a procedure is 
without problems if one sticks to the collinear dynamics. In 
Sec. |III]below, however, we include primordial transverse mo- 
menta of the partons, i.e. we have to deviate from p, = x,P,. 
The natural choice would be to define x, via one nucleon mo- 
mentum component (the large one). This is exactly what we 
have done here for the collinear case. In the naive parton 
model x, and x,, i.e. pf and pj, are introduced as independent 
variables which are then constrained by the kinematical and 
onshell conditions (T30b-(f33T)- However in the Bjorken limit 
(M, S — > oa,M 2 /S = const) the parton momenta should be- 
have like 01 



P \ = 0(M), p\ = 0(1 1 M), 
pl = (KM), p- 2 = 0(1 1 M). 



(47) 
(48) 



For x, + this power counting is not fulfilled, cf. Eqs. d27l ) 
and d29l . Hence this solution corresponds to non-factorizing 
power suppressed corrections. Therefore in the naive parton 
model one falls into a trap by picking up this additional un- 
physical solution. The same happens for the more compli- 
cated case including primordial transverse momenta. 
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It is worth pointing out the connection between the two 
types of parton models (standard vs. naive) and QCD. There, 
e.g., the DY cross section formula emerges from factorisa- 
tion. It turns out that in the Bjorken limit a PDF depends on 
one variable only 13211 . which is encoded in x,. In the final 
formula the energy-momentum relation, e.g., for the DY pro- 
cess takes the form 6(q - x\P\ - X2P2) which suggests the 
interpretation of as the parton four-momentum. Thus the 
standard (collinear) parton model emerges from QCD and not 
the naive one. 

Including in addition primordial transverse momenta one 



dcr na j ve 
dM 2 dx F 



The unphysical second solution for the momentum fractions 
is represented by 

S(x 1 )6(x 2 ) xJiixuM 2 ) X2MX2, M 2 ) (50) 

in the last expression. Its contribution does not vanish since 
one obtains for large enough M 2 ll33ll 

lim (x/(x,M 2 ))>0. (51) 



Note that in the last expression we have recovered the standard 
parton-model result Eq. dZ5b . 

The main reason to present this naive approach in detail will 
become clear in the next section where we lift the simplifica- 



has to model the distributions of these momenta. However 
there is a constraint the chosen model has to obey: in the 
Bjorken limit one should come back to the standard parton 
model and not to the naive one since only the former emerges 
from QCD. 

In the following we will point out how to modify the naive 
parton model such that one ends up with the standard parton 
model. This procedure will then be generalized to the case 
where primordial transverse momenta of the partons are in- 
cluded. In the naive parton model the hadronic cross section 
reads: 



I 

We now introduce a notation which we will keep through- 
out this paper. Whenever we explicitly disregard unphysical 
solutions of the type of Eqs. (f39l> - (|43l> under an integral we 
denote this integral by f. Thus 



(52) 



(53) 



tion of a collinear movement of the partons with the nucleons. 



= JT dJCi JT dx 2 J] qfMxi , M 2 ) Mx 2 , M 2 ) ^ T s(M 2 -(p 1 + p 2 ) 2 ) 6 {x F - 

P' a C X a V JZ*f a A 2 47ra 2 2(q z ) m . ix x l x 2 (S(x l - x l+ )5(x 2 - x 2+ ) + 5(xi)6(x 2 )) 

= d*i dx 2 X ^Mxi,M 2 )Mx 2 ,M 2 )— — — — (49) 

Jo Jo • 9M 2 S i/Z \(x\X2 - Xix 2 ){x\ + x 2 + X\ + x 2 )\ 



dcr naive f 1 f' . V 2,, .j-z^ j-, , .2 4na 2 2(q z ) max x l x 2 (6(x l - x l+ )5(x 2 - x 2+ ) + 6(x l )6(x 2 )) 

— = dxi dx 2 > qf fi{x\, M ) fi(x 2 , M ) T — — 

dM z dx F Jo 'Jo 2 Zj t/ ' /,v ;M 9M 2 S^ 2 \(x l X2-x l X2)(xi+x 2 +x l + ~X2)\ 

ras 

do- f 1 , f 1 y" 1 2r/ ,a ,, » ,2 ,4na 2 2(q z ) max xix 2 {5{xi - x u )6(x 2 - x 2+ ) + 6(xi)6(x 2 )) 

— — = + d Xl + dx 2 X qifi(xi,M 2 )fi(x2,M 2 )— 2 ^nV ( * ~ * 1 ~ m 

dM z dx F Jo Jo 9M Z S ilz \(xiX2 - xix 2 )(xi + x 2 + X\ + x 2 )\ 

- X' * X' ^mm,m^ 2( ^^ ) -;^-- ) 

Z2r, ,,2s r, „ ,2s Atllal ("?z)max 
q z fi(x u ,M z ) f-j(x 2t ,M ) t T, • 



max 
. +*2 + ) 
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III. FULL KINEMATICS 

The Bjorken limit and the corresponding infinite- 
momentum frame in which the standard parton model is well 
defined and derived from leading-order pQCD is an idealiza- 
tion of real experiments. There the nucleons will always move 
with some finite momentum and thus the partons inside the 
nucleons can have nonvanishing momentum components per- 
pendicular to the beam line. The factorisation into hard (sub- 
process) and soft (PDFs) physics is proven in the collinear 
case at least for leading twist (expansion in 1 /M) in 113411 and 
in the transverse case at least for partons with low transverse 
momentum in ll3~5ll . 



A. Transverse-momentum distributions 

For the calculation of the hadronic cross sections we 
will need transverse-momentum dependent parton distribution 
functions. We denote these by f t . They are functions of the 
light-cone momentum fraction x/, the transverse momentum 
Pi ± and the hard scale of the subprocess q 2 . The general form 
of these functions however is unknown. Known rather well 
are the longitudinal PDFs. Since data of DY pair production 
are compatible with a Gaussian form of the /^-spectrum up 
to a certain pj IL36L l37ll . we assume factorisation of the lon- 
gitudinal and the transverse part of fi and make the following 
common ansatz ll20il38ll39ll 

Mxi,Pi ± ,q 2 ) = /KW) ' MPiJ ■ (54) 

Here fi are the usual longitudinal PDFs and for fi ± we choose 
a Gaussian form, 

^ ) = i^ eXP (-w)- (55) 
I 



In this approach the transverse momentum (pj = \q ± \) of the 
DY pair is accessible, since the annihilating quark and anti- 
quark can have finite initial transverse momenta. Note that in 
the calculations of J20ll the partonic DY cross section dcr was 
taken in the collinear limit. 

I 



The width parameter D is connected to the average squared 
transverse momentum via 

{(Pl) 2 } = Jdjf i Jjf i J 2 f i APiJ=4D 2 (56) 

and it has to be fitted to the available data. 



B. Cross section 

Now we calculate the hadronic cross section dcr taking into 
account the full kinematics. Since it is necessary to remove 
the unphysical solutions for the light-cone momentum frac- 
tions x\ and X2 which correspond to the ones found in Sec. 
Ill B I for the collinear case, the calculation has to be conducted 
such that it is possible to disentangle the physical and the 
unphysical solution. First we will discuss in Sec. IIII B 1 1 a 
straightforward calculation which however does not obey this 
requirement (the details of this calculation can be found in the 
Appendix). If one does not remove the unphysical solutions 
one produces unphysical results. This reveals a pitfall which 
the unawareness of this problem can create II27I429I1 . In Sec. 
IIII B 21 we will show how to properly remove the unphysical 
solutions as we did for the collinear case at the end of Sec. 
HTBl 

In the transverse-momentum dependent approach the 
leading-order Drell-Yan total differential cross section reads 
10 



(57) 

I 

1. Naive calculation 

In the naive approach the partonic triple-differential cross 
section reads: 



- = f * K - 1 2 ) ^ + P2 - *) * U - ^As {p 2 t - (4,) 2 ) 

= ^^K-to + ^)V(^-l^)^(^-^+Aj 2 ) ■ (58) 
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Inserting Eq. d58b in Eq. $57\ yields a multiple integral for the 
triple-differential cross section: 



dcr„ 



dM 2 d.v, vi 



= | dxi I dx 2 I dp I± I dp 2± F(x u p l± , x 2 , p 2± , M 2 ) 
F apj Jo Jo J J 

xS(M 2 -( Pl + P2 ) 2 ) 6 (x F - (Pl - Pl)z ) 6 (p 2 T - (p 1± + frj 2 ) . (59) 



All pieces which do not contain ^-functions are collected in 
F{...). The straightforward, but naive calculation of d59l ) was 
performed in |27T - I29tl . The details of this calculation can be 
found in the Appendix. Here we just want to point out the 
problems arising from this approach: the naive calculation 
with the full kinematics incorporates unphysical solutions for 
the momentum fractions x, which correspond to the unphysi- 
cal solutions of the collinear case in Eqs. (|39l-(l43ll. However 
in the collinear kinematics it is quite clear that these solutions 
for Xi cannot be the physically interesting ones, since they are 
just Xi = and the PDFs are divergent for small x as one can 
conclude from (IBTT ). In the case of full kinematics the situa- 
tion is similar, however due to the introduction of transverse 
quark momentum distributions the momentum fractions Xj are 
smeared out around their collinear values. Nonetheless the 
unphysical solutions are still very close to zero and one picks 



up very large contributions of the diverging PDFs at such low 
x. This leads to a large enhancement of the cross section in the 
full kinematical approach and the data are now overestimated. 
The effect can be seen in I27H29I1 and also in Sec. II VI where 
we compare this naive approach and the correct calculation of 
the next section. In addition, in Sec. ||V]it will be shown that 
the M-dependence of the cross section is not reproduced in 
the naive approach. 

2. Correct calculation 

The analogue to Eq. ( |59l in the correct approach is (for the 
notation J- see Sec. IHBl l 



— — T = -f dxi -f dx 2 f dpi ± f dp 2± F(xi,pi ± ,x 2 ,p 2± ,M 2 ) 
F dpj Jo Jo J J 

x6(M 2 -( Pi +P2 ) 2 ) sL - ( Il±Elk) 5(p 2 T - (P 1± + p 2 J 2 ) 



The (5-functions in Eq. (f60b must be worked out in a way that 
allows to discern physical and unphysical solutions for the 
momentum fractions x, in order to perform the ^--integrations. 
For this aim it is useful to rewrite the parton momenta in terms 
of different variables: 



Adding and subtracting Eqs. (l63l l and d64b yields 



k 2 = --M 2 , (65) 
4 

k ■ q = . (66) 



q = pi + pi , 

k = ^(p2 - pi) ■ 



(61) 
(62) 



Inverting the last two equations, we can use the on-shell con- 
ditions for the partons to get 



and 



= p\ = \-q-k \= -q L -k-q + k L (63) 



= p\ = ( i q + k\ = X -q 2 +k- q + k 2 . (64) 



Solving Eq. ( I65l l for k + yields 



h+ _ k 2 ± -\M 2 



(67) 



Inserting this result into Eq. d66l > gives an equation quadratic 
in k~: 

= k + q~ + kT q + — 2k ± • q* ± 
t 2 - iM 2 

= — A - q~ + k-q + - 2k± ■ q ± (68) 

k 

= (r ) 2 q + - 2k ± ■ 4JC + - \m 2 \ q- . (69) 
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The solutions are 



k, ■ q. 



k± ■ q± 

q + 



(70) 



Inserting (|70l > into d68t gives the solutions for /: 
k± ■ q± 



q 



k± ■ qi 

q + 



(71) 

Rewriting now Eqs. (|26T > and (l28l i in terms of g and k we obtain 
the solutions for the parton momentum fractions: 



(*i)± 



9" - (n± 



and 



(* 2 ), = -^z = ^(^ + + (fe + ), 



(72) 



(73) 



Since there are two solutions for k~ and k + , respectively, we 
also get two solutions for x\, x 2 . To determine which set of 
x\,X2 and thus k + ,k~ has to be chosen we take the limit of 
zero parton transverse momentum. In this way one can make 
the connection to the collinear case (then q 2 — > q + q~ - M 2 ): 



(74) 
(75) 



Inserting expressions ( |74"1 i and d75l l into d72b and ( l73l yields 
two solutions for the momentum fractions, just as in the 
collinear case in Sec. Ill Bl 



and 



(*i)i 



-1° 



- 

V5 W + 



(76) 



(77) 



The lower solutions correspond to the standard parton model 
Eqs. (05), (ED, since XiX Z = ^andx 2 -*i = % = x F ^f^. 
The upper solution then corresponds to the unphysical case 
x\ = X2 - and jci # + xi, see Eqs. d39b-d43l. 

This is the crucial point: to receive physically meaningful 
results from Eq. ( f60b one has to discard these upper solutions 
just as one does in the collinear case in Sec. Ill Bl This requires 
that the integrals in Eq. d60l ) are evaluated in the correct order, 
otherwise one cannot disentangle the two different solutions 
for X] and xi- We will now present a calculation which re- 
spects this requirement. In Sec. [IV] we will show that the 
quantitative difference between this calculation and the calcu- 
lation from Sec. IHIB H is huge. 

We begin by introducing several integrals over (5-functions 
in Eq. d60l >. In this way we will transform the integration vari- 
ables to the above chosen q and k ± : 



dM 2 dx F dp 



j- dxi^ dx 2 J dpi ± J dp* 2x J dq ± J dk ± J dq + J dq F(x u pi ± ,x 2 ,p 2± ,M 2 ) 
X 6(q + - {p\ + p + 2 )) 5{q- - (p\ + pi)) # 2 \q* ± - {p lx + 2 J) 5 {2 \h - \ (j} 1± - pi j) 



First we perform 



x S (M 2 - ( Pi + P2 ) 2 ) 6 [x P - (Pl - P2)z )d (p 2 T - (p l± + p 2 J 2 ) . (78) 
j dp l± J dp 2± 5 (2) (q* x - (p l± + p 2 J) (5 (2) fe x - ~ (p 2± - p 1± )\ = 1 . (79) 



I 

Now we calculate the integral According to Eqs. (I70l>-(l73l the <5-functions in the last expres- 

sion have two possible solutions for each p~ and pi,. However 
■f 1 d -f 1 A x< + ' + +Vl x< ~ ' ~ as explained above we now have to explicitly remove the un- 

Jo Jo 



dx 2 6(q (p 1 + p 2 )) 5{q (p l + p 2 )) . physical solutions (xi)+ and {xi)- , which are the ones corre- 



(80) sponding to the upper sign in Eqs. ( 1701 and ( TTTb : 
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f d*, f 
Jo Jo 

Jo Jo 



djc 2 c% + - OJ" + p 2 )) - (Pi + P2)) 



dx2 6 



<? — x 2 VS 



*1 



<7 -*i 



= J dxi J dx 2 5U ^--fe)+V5 



-(jci)_Vs - 



(Kt+jx) n 



2\ 



(3^ - kj\\4x + 



5 - 



( Xl ) 2 _(x 2 )lS 

Using dg+di^r = 2d^od^ ; we can evaluate some of the remaining integrals of Eq. ( 1781 with the help of the (5-functions: 

j dq + dq-d4 ± 6 (M 2 - q 1 ) 5 [x F - J^—) * {pt ~ (4±f) 



= 2 J dq d4 ± dq z 6(M 2 + (4 ± f + q\ - qfj 5\x F - J^—j - (^) 2 ) 



* (<7;)n 



with E = -,Jm 2 + pj + xpiq^Bs^ ■ Collecting the pieces, what remains of Eq. d78l is 



da = 

' dx F dp 2 T J 



dM 2 dx f 

(x\)-,p\ x , (x 2 )+ and p\ x are now fixed: 



dk ± 



(piJHpiJ 2 



s - 



(*l)_(* 2 )^ 



(Xl)- = 


1 


q 
2 


(xi)+ = 


1 

vf 


2 


fax = 


1 . 




fax = 


+ £x 




with 



/max 
max 



q + = E + x F (q z ) B 
q~ = E - x F (q z ) B 

\q±\ = Pt , 
k±- q± = \lc±\pr cos(0 x ) . 

I 



(81) 



(82) 



(83) 



(84) 

(85) 

(86) 
(87) 



(88) 
(89) 
(90) 
(91) 



|£j_lmax is fixed by the condition that (xi)_ and (jc 2 )+ must be 
real numbers: 



(k ± ) 2 < 



We have convinced ourselves that this condition also guaran- 
tees that < (x\)-, (x 2 ) + < 1. Finally we arrive at the follow- 
ing expression: 



M 2 + p\{\ -cos 2 (^ ± )) 



(^)L X (92) 
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AM 2 dx F dp\ Jo X Jo 



2^ £ 



(/?lJ 2 (p 2 J 2 



5 - 



( Xl y_( X2 )is 



F(( Xl ).,p K ,(x2) + ,P2„M 2 ) (93) 



with 



Ana 2 



i 

and with // defined in Eq. 



(94) 



IV. RESULTS 

In this section we present our quantitative results and com- 
pare the naive approach of Sec. IHIB II and the correct ap- 
proach of Sec. IHIB 21 The data are from the NuSea Collab- 
oration (E866) |36j|37ll and from FNAL-E439 H. For the 
collinear PDFs we used the GRV98 LO parametrization ll33tl 
available through CERN's PDFLIB version 8.04 J4lll . 



A. E866 - p T -spectra 

Experiment E866 measured continuum dimuon production 
in pp collisions at S ~ 1500 GeV 2 . The triple-differential 
cross section as given by the E866 collaboration is 



dcr 



2E dcr 



d 3 P ny[S dx F dp 2 T 



(95) 



where an average over the azimuthal angle has been taken. 
The data are given in several bins of M, x F and p F and for ev- 
ery datapoint the average values (M), (x F ) and (p?) are given. 
Since our schemes provide Eqs. ( IA8t and d93l we calculate the 
quantity of Eq. d95l l for every datapoint using these averaged 
values and then perform a simple average in every M 2 -bin: 



2E dcr 
7T VS dx F dp 



2E 
2E 



f 



dcr 



M 2 -bin dM 2 dx F dp 2 

dcr 



-dM l 



„ V5 dM 2 dx F dpt 



■((M),(x F ),(Pt)) 



where 



= a/<M> 2 + (p T ) 2 + {x F ) 2 {(q z ) a 



(96) 



(97) 



with M max (M m i„) the upper (lower) 



and AM 2 = M 2 - M 2 . 

max mi 

limit of the bin. 

We plot the results for the two different approaches in differ- 
ent M-bins in Fig. |2] Everywhere a value of D — 0.5 GeV for 
the transverse momentum dispersion was chosen. The solid 
lines represent the correct approach. The shape of the spectra 
is described rather well which is due to the choice of the pa- 
rameter D. However the absolute size is still underestimated 
and a factor K m 1.75 - 2.0 would be necessary to repro- 
duce the height of the data. The naive approach is plotted 



with dashes. As already mentioned in Sec. IIII B 1 1 the calcu- 
lated cross section overestimates the data significantly. This 
can also be seen in B27I429I1 . We note that the discrepancy be- 
tween both approaches is about 1 order of magnitude and it 
becomes worse in the higher mass bin. This already indicates 
a wrong M dependence of the naive approach. 
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FIG. 2. pr-spectrum obtained from the naive and the correct ap- 
proach with D = 0.5 GeV. Data are from E866 binned with 
4.2 GeV < M < 5.2 GeV and 7.2 GeV < M < 8.7 GeV, 
-0.05 < Xf < 0.15. Only statistical errors are shown. 
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B. E866 - M-spectrum 

The double-differential cross section is given by the E866 
collaboration as 



M 



Ao- 



AM Ax F 



(98) 



Again the data are given in several bins of M and Xf and for 
every datapoint the average values (M) and (xp) are provided. 
Once more we start with Eqs. dA81 > and ( l93l and calculate the 
quantity of Eq. ( |98l by integrating over p 2 T for every datapoint 
using these averaged values: 



da 



AMAx F 



(My 



Ap\ 



Act 



AMAx F Ap 2 



= {M)- 







Ap\ 2 (M) 



Act 



AM 2 Ax F Ap 



«M>,<x F » . 



(99) 



The maximal possible pj is determined by the kinematics. 



P x + p 2 = q + X 
(Pi +Pi- qf = X 2 = M\ 
> S + M 2 - M\ = 2 (Pi + Pi) q 
= 2 y[SE 



M 2 + (Pt)L x 



= 2 Vs" ^M 2 +p 2 T +q 2 
(S + M 2 - M 2 ) 2 



+ q 2 z =E 2 = 



(Pr) ma , 



45 

(S + M 2 - M 2 ) 2 
4S 



(100) 
(101) 

(102) 
(103) 



M L 



(104) 



M 2 is the minimal invariant mass of the undetected remnants. 
We choose a value of M« = 1.1 GeV. Note that at cm. ener- 
gies of y/S ~ 27.4 GeV (E439) and V5 * 38.8 GeV (E866) 
we are not really sensitive to this value if it stays at or below 
a few GeV. 

The results are plotted in Fig. [3] Again we use D = 0.5 
GeV. The solid line represents the correct approach, the long 
dashed line the naive one. For comparison the result of the 
standard (collinear) parton model is plotted with the short 
dashed line. Here the discrepancy between the naive and the 
correct approach is fully visible, since neither the slope nor 
the size of the M-spectrum is reproduced in the naive ap- 
proach. Instead it gives almost a constant distribution. (Note 
here that this dataset is not shown or compared to calculations 
in 1I27I - I29I1 .) The correct approach however describes the slope 
well and again a factor K x 1.75 - 2.0 is necessary to reach 
the absolute height of the data, as expected from the triple- 
differential results in the last section. Note that the result of 
the correct approach and the standard (collinear) parton model 
coincide. 



N " 100 
> 

<D 

CD 

JO 

r 10 





"O 



0.1 




correct 
naive 
parton model 
E866 



7 8 9 
M [GeV] 



10 11 



FIG. 3. M-spectrum obtained from naive and correct approach with 
D = 0.5 and from the standard parton model. Note that the results 
from the latter two approaches are basically on top of each other. 
Data are from E866 binned with -0.05 < x F < 0.05. Only statistical 
errors shown. 



C. E439 - M-spectrum 

Experiment E439 measured dimuon production in pW col- 
lisions at S ~ 750 GeV 2 . The double differential cross section 



da 



AMAx' 



has been given with 



x F 



1 



M 2 
S 



(105) 



(106) 



at a fixed x^. = 0.1. 

As before we begin with Eqs. ( IA8I ) and d93l and calculate 
the quantity Eq. ( 1105b by integrating over p 2 T and performing 
a simple transformation from Xp to x' F : 
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-= dpi 7 = ApllM 1-— — 7 \M,x F = x' F 1-— . (107) 

£ Jo lT dMdx' F dp 2 T Jo ^ r \ S dM 2 dx F dp\\ F \ S 



We plot the results in Fig. |H the solid line represents the 
correct approach, the long dashed line the naive one. With the 
same parameter D = 0.5 GeV as for the E866 case we find 
the same discrepancy between the two approaches. Again the 
correct approach reproduces the slope well and a factor K » 
1 .6 is required to fit the data, while the naive approach fails to 
describe the slope and absolute size of the cross section. Once 
more the result of the correct approach agrees well with the 
result of the standard parton model (short dashed). 

Here we note the following: in II27I429I1 the same data of ex- 
periment E439 are compared to calculations, however only to 
an approach including both initial quark transverse momen- 
tum and quark mass distributions. There it is found that the 
data can be described well without a K-factor. There is no 
comparison of E439 data with a transverse momentum de- 
pendent calculation with onshell quarks (i.e. what we call the 
naive approach) in 1I27I - I29I1 . We acknowledge that the in- 
troduction of quark mass distributions lowers the cross sec- 
tion which somewhat compensates for the enhancement in 
the naive approach. Nonetheless we want to point out that 
even with additional smearing from the quark mass distribu- 
tions the naive approach will always lead to the wrong ./In- 
dependence of the cross section. The reason is simply that the 
PDFs are probed in two areas: around the standard collinear 
parton model x, cf. Eq. ( f25l l, and in a region very close to 
x = where the PDFs behave very differently with M and 
give much larger contributions than for the physical x, since 
the PDFs diverge rapidly for x — > 0. Thus we conclude that 
the agreement of the full calculation with the E439 data in 
127-29] must be erroneous. 




4.5 5 5.5 6 6.5 7 7.5 8 
M [GeV] 

FIG. 4. M-spectrum obtained from naive and correct approach with 
D = 0.5 and from the standard parton model. Note that the results 
from the latter two approaches are basically on top of each other. 
Data are from E439 with x' F = 0.1. 
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V. CONCLUSIONS 

In this paper we reexamined a phenomenological model of 
Drell-Yan pair production II27I429I1 . This model tried to de- 
scribe the DY process in a parton model scheme which takes 
into account the full transverse parton kinematics in the hard 
subprocess. The aim was to reproduce the transverse momen- 
tum spectra of the DY pairs and in addition to reproduce the 
absolute size of the cross sections by introducing mass distri- 
butions of the partons. 

We have shown that already in the first step of introducing 
full transverse kinematics important constraints were not con- 
sidered. Unphysical solutions emerging in a (too) naive par- 
ton model contaminate the results. We have derived these con- 
straints in the usual collinear approach and then made the con- 
nection to the more general case of full kinematics. It turned 
out that unawareness of these constraints can lead to a dras- 
tically different result of the calculations: while in j27l - l29ll 
the inclusion of the transverse kinematics in the subprocess 
leads to an overshoot of the cross section, our corrected ap- 
proach shows no such behavior and instead nicely reproduces 
the standard parton model prediction for the invariant mass 
spectra and the low transverse momentum spectrum, however 
only up to a K factor. Additionally we find that the naive 
approach taken in ll27l - l29ll produces a wrong M-dependence 
of the cross section. This is a crucial point since already the 
standard parton model reproduces the right slope of the M- 
spectrum. Therefore we conclude that the findings in I27M29I1 
that allow for a K factor free description of DY pair produc- 
tion are unwarranted. 
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Appendix A: Naive calculation of the hadronic cross section 
with full kinematics 

Rewriting the ^-functions of Eq. d59l in terms of the inte- 
gration variables yields 
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S(M 2 - (pi + pi) 1 ) di^Xp 
■■S{m 2 -pIp2 -PiPt + 2pi x -p 2x ) s(x F 



(P\ + P2)z 
(<7z)max 



s(pt~(Pi x +P2J 2 ) 

P + l ~ Pi + Pi ~ Pl\ ,( 2 fit .2 ,* * * \ 



Xf - 



J M 2_ iP±L (^L _ p - pt + p 2 _ ^ja _ (/?2j2 ) tf 
\ ^1 ^2 / 

Xc5(pr-(/?ix) 2 -(P2j 2 -2/?i ± -^2 ± ) 

J M > - _ Xl V5*2 Vs +P 2 - (A J 2 - (P2J 2 ) 6 

\ x x yfS x 2 Vs / 
x 6 (p\ - (p lx f - (p 2 J 2 - 2p lx ■ p 2x ) . 



x F ■ 



2(q z ) n 



(Al) 



Note here that if one puts all transverse momenta in Eq. ( IAU to zero the collinear relations (I32I33I I are recovered. Note also that 
the unphysical parts of Eqs. ( f32l > and ( 1331 are included here since 



Xl = 



x 2 



p\ ; (plj 2 

Vs nVs" 

/>2 (P2j 2 



Now using the first two 5-functions of Eq. (IAU we can obtain solutions for the squared transverse momenta: 



(A2) 
(A3) 



SIM- 



2 (PO 2 (Plj 2 

xi V5 x 2 VS 
X1X2S 



xi Vsx 2 Vs + /4 - (pij 2 - (/? 2 J 2 ) ^ 



xf • 



2(<7z)„ 



= (^)max "~77 (5 ((Plj 2 -(E-X 2 yfS+ (q z )m^X F ) X X Vs) (G?2 ± ) 2 - (j5 - Jtl VS - (<7 Z )maxXf) X 2 Vs) 



with the energy of the virtual photon 



E = ^M 2 +p 2 T + (q z ) 2 mm x 



(A4) 



(A5) 



We note that exactly at this point the physical and the unphysical solutions for the momentum fractions x, have been mixed up 
by rewriting the <S-functions, since the unphysical solutions of Eqs. ( IA2b and dA3b have entered. 
Transforming the transverse momentum integrals 



dpi ± J dp 2± = J d0 lx J dfe J 2 d(/?u)2 J 2 d(/?2j 



(A6) 



we can rewrite the entire expression d59t in the following form: 



do~naive 



dM 2 dx F dp 



r*\ r*\ r*7T r*2.7l r*oa ^ *-»co 

- = dx, dx 2 d<p ix d4, 2± -d(PiJ 2 -d(p 2± 
T Jo Jo J-7T Jo Jo z Jo L 



F(xi,pi ± ,x 2 ,p2 ± ,q 2 ) 



, (q Z )m!xxXlX2S 

E ^(P0 2 (P2j 2 ~ (p\ ~ (Plj 2 ~ (P2J 2 ) 2 

: 6 ((pi J 2 -(E-x 2 ^S+ (q z ) m ^x F ) xi Vs) 5 ((p 2 J 2 -(E-Xi^/S- (q z ) m , x x F ) x 2 Vs) 
x 5 <t>\ x - arccos 



P 2 T ~ (PiJ 2 ~ (P2J 2 



(Al) 
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Now all four integrations concerning the partons' transverse momenta can be carried out, leaving a two-dimensional integral 
which must be calculated numerically: 

do- naive f r 7r(q z ) mlix x 1 x 2 S F(x u p l± ,x 2 , p 2± ,q 2 ) 

r- = I d^i I ax2 — . (A8) 

dM*dx F dp 2 T J J E I 2 2 , 2 2 ,2 

y4( Pli ) {P2 ± ) -(p 2 T -{pi ± ) -{?2 X ) ) 

(p\^ and (p2 ± ) are given by the 5-functions in Eq. ( lA4b and the integration boundaries of x\ and x 2 have to be chosen such 
that the requirements 

CpiJ 2 > , (A9) 
CpiJ 2 > , (A10) 
MpOHhj 2 ~ {pi ~ CpO 2 ~ CpiJ 2 ) 2 > (All) 



are fulfilled. One finds 



E x F (q z ) m3X 

< xi < — (A12) 



and 



*1 



VS (M 2 - p\) + p\(E- x F (q z )™x) ~ 2Mp T ^ V?(£-M? Z )max-^I VS") 

VS(£ - jc f fe) max ) 2 

< JC 2 < 

xi VS" (M 2 - p 2 ) + p\ (E - x F (q z ) max ) + 2Mp T yfxi ^(E - x F (q z ) max - x x y/S) 

y/S(E - x F {q z ) max ) 2 



(A13) 
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